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Abstract 



In this paper the Levy-Leblond procedure for linearizing the 
Schrodinger equation to obtain the Pauli equation for one particle 
is generalized to obtain an iV-particle equation with spin. This is 
achieved by using the more universal matrix factorization, GG = 
\G\I = {—K) l I. Here the square matrix G is linear in the total energy 
E and all momenta, G is the matrix adjoint of G, I is the identity ma- 
trix, \G\ is the determinant of G, I is a positive integer and K = H — E 
is Lanczos' extended Hamiltonian where H is the classical Hamilto- 
nian of the electro-mechanical system. K is identically zero for all 
such systems, so that matrix G is singular. As a consequence there 
always exists a vector function 9 with the property G9 = 0. This fac- 
torization to obtain the matrix G and vector function 9 is illustrated 
first for a one-dimensional particle in a simple potential well. 

The quantization of the matrix relation, G9 = 0, is performed, 
using the Planck-de Broglie laws and a Fourier transform to obtain 
the generalized wave equation, G(p, q)tp(q) = 0, where q and p are, re- 
spectively, the generalized position and momentum-operator vectors of 
(n + l)-dimensions. This quantization procedure, when applied to the 
above one-dimensional example, yields immediately a pair of coupled 
operator equations which become, after eliminating one component of 
ip, the classical second-order partial-diffferential Schrodinger equation 
of a single one-dimensional particle in a force field. This same tech- 
nique, when applied to the classical nonrelativistic Hamiltonian for iV 
interacting particles in an electromagnetic field, is shown to yield for 
N = 1 the Pauli wave equation with spin and its generalization to N 
particles. 

Finally this nonrelativistic generalization of the Pauli equation is 
used to treat the simple Zeeman effect of a hydrogen-like atom as a 
two-particle problem with spin. This analysis of the Zeeman effect in a 
weak magnetic field exhibits the usual two-fold splitting of the energy 
levels, obtained first by Pauli using the Pauli equation for N = 1, 
except for a slight change in the particle mass needed in the Larmor 
frequency. This "increased" mass, given by rriL = m\mil {vn,2 — mi), 
agrees with the g-factor correction for nuclear motion obtained in 1952 
by W.H. Lamb using the relativistic, one-particle, Dirac equation with 
potentials. Also it is consistent with the known fact that positronium 
has no weak-field Zeeman effect. 
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1 Introduction: Extended Hamiltonian and 
its Factorization 



The classical Hamiltonian 

H = H(q u - ■ ■ ,q n ;pu- ■ ■ ,p n \t), (1) 
which is also called the "total energy," is related to the Lagrangian, 

L = L(q u ■■■ ,q n ;qi, ■ ■ • ,g„;t), (2) 

by the relation, 

n 

H = ^2p k q k -L, (3) 

k=l 

where the conjugate variables q k and pk for (A; = 1, 2, • • • , n) are, respectively, 
the generalized coordinates and momenta of the particles of a mechanical 
system in a force field. The momenta p k are obtained from the Lagrangian 
by means of the formulae, 

p k = dL/dq kl (4) 

for k = 1, . . . ,m in terms of the time derivatives q k of their conjugate co- 
ordinates q k . The Hamiltonian H in ([!]), as a function of the momenta p k , 
is obtained from (j4]) by solving the equations (jl]) for the q k , followed by a 
substitution into ([3]), e.g., see [1, Chapt. 6] and [2, Chapt. 8]. 

In classical mechanics a path minimization of the action integral, 

A = I AnAxr-- ,q n ;t)dt, (5) 

yields Lagrange's equations, a system of n ordinary differential equations for 
the "path" (<?i(i), 92^), • ' ' > Qn(t)) of the system point as a function of time 
t. However, sometimes it is advantageous for this minimization process to 
elevate time t to the position of the other coordinates q k . This means that 
time is taken to be some monotonically increasing differentiable function, 

t = t(a), (6) 

of some real variable parameter a. 



3 



In terms of the parameter a in ([2D the action A in (jSJ) becomes, by an 
evident change of variables, the integral, 

A = I L(Qir--,Qn;q'i/t'r--,q'Jt';t)t'da, (7) 
Jti 

where q±, - ■ ■ , q n , t are functions of a, q' k = dqk/da, and t' = dt/da. In (J7J) the 
original Lagrangian in the integrand of ([5]) clearly becomes the new modified 
Lagrangian, 



L{qi,--- ,qn,qn+r,q'i,--- ,q' n +i) 

L(qx,-~ , q n ; q'i/t', ■ ■ ■ ,q'Jt';t)t', 



(8) 



where q n+ i = t. Formula (0J is applied next to L in (|SJ) to obtain the 
momentum associated with time t, namely, 

Pt = dl/dt' 

= L + (^dL/dq k dq k /dt^jt' 

= L ~ {^dL/dq.q'Jt'^t' 

' n \ 

^Pkqk-L . (9) 



K k=l 



Hence the momentum p t is the negative of the total energy defined in ([3]), 
e.g., see [1, Sec. 5.8] and [2, Sec. 8-4]. Note that if the Lagrangian L is not 
an explicit function of time t, then the system is conservative. In this case 
H = E, where E is a constant. Thus the momentum p t in (Q is given by 

Pt = ~H = -E. (10) 

If the system is nonconservative, it is convenient to still let E denote the 
total energy. Thus in all cases the relation in fllOl) . between the momentum 
p t and the total energy E, must hold. 

Now let the new momentum p t in ([9]) be denoted by 

Pn+l = Pt, (11) 
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the momentum conjugate to the time coordinate q n+1 = t. Also use ([2D to 
obtain the Hamiltonian, 

n+1 

H = Y,Pkq' k -U (12) 

fe=i 

in terms of the modified Lagrangian, 

L = Lq' n+1 , (13) 

the derivatives q' 1: q' 2 , ■ • • , q' n+l with respect to parameter a, and the momenta 
Pi,P2, ■ ■ ■ ,Pn+i, where p n +i is defined in f fTTT) . A substitution of L, given in 
( TIB"]) , into (fT2"l) . using the Lagrangian L in Q, yields finally the formula 

n+1 



# = ^Pkq'k - 5^ p ^ fe _ H ) q 'n 



k=l 

= ( Pt + H)t' = Kt', (14) 

for the modified Hamiltonian in the parametric system as a function of pt 
in fflOl) and the original Hamiltonian H in ([1]). By ([TO]) and the assumption 
that t'(a) > for all a the expression in ffl4l . namely, 

K=(p t + H) = -E + H (15) 

is always the constant, zero, for any electro-mechanical system and is called 
by Lanczos in [1, pg. 190] the extended Hamiltonian function. Since K 
does not depend explicitly on the independent variable a, every system can 
be considered to be a conservative system of the extended phase space of 
2(n+ l)-dimensions. Lanczos further points out that any "fluid particle" of 
this phase space is "steady" and permanently on the surface, 

K = const. = 0, (16) 

of the phase space. In either a classical or quantum mechanical system the 
identity in (flBT) . satisfied by the extended Hamiltonian in ffTBT) . constitutes 
the "fundamental information" of such systems. 

In quantum mechanics the relation in fflBI) is used next in a manner similar 
to the matrix factorization used by Dirac [3, pg. 255] to obtain his famous 
relativistic wave equation. To accomplish this we note first by ffl5|) that 
the momentum, p t = —E, occurs linearly in K whereas all other momenta 
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Pi) • ' " iPn generally occur quadratically. To put the momenta pi, • • • ,p n and 
Pt on the same footing it is desirable, following Dirac [3] and Levy-Leblond 
[4] later in 1963, in any matrix factorization of K that all momenta occur 
linearly. 

In order to achieve a matrix factorization of K in such a manner that the 
momenta occur linearly use is made next of the well-known matrix identity, 

GG=\G\I, (17) 

where \G\is the determinant of some square matrix G of complex numbers, G 
is the matrix adjoint of G, sometimes denoted by Adj(G), and finally I is the 
identity matrix. It is well known that the adjoint of G, namely G = Adj(G), 
equals the matrix transpose of the matrix of cofactors of G, e.g. see [5, Sec. 
1.6]. Thus the universal factorization identity in ffTTj) is closely related to 
Cramer's rule for solving a system of linear equations. 

Here, using ffT7|) . a matrix G is said to linearly factor the extended Hamil- 
tonian K if the following two conditions hold: 

i) . Matrix G is linear in the momenta p±, ■ ■ ■ ,p n and pt- 

ii) . \G\ = (-K) where / is a positive integer. 

To illustrate the linear factorization in i) and ii) an example of the ex- 
tended Hamiltonian for a one- dimensional particle in a force-field is given 
next. 



Example 1 Let a one-dimensional particle of mass m be in a time-varying 
force field of potential V(x,t). Then the Hamiltonian is 



H(x, t; p) = p 2 /2m + V(x, t), 



(18) 



where p is the momentum of the particle. Thus by (j!5p and f)16p the negative 
of the extended Hamiltonian for this system is given by 



-K{x,t;p, -E) 



E - H(x,t;p) 
E-V - p 2 /2m 
0. 



If one lets G in ([17]) be the 2x2 matrix 
G = 



E-V p/V2m 
pj \Jlm 1 



(19) 



(20) 
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then its determinant is computed easily to be 



\G\ = -K(x,t;p, -E) = 



(21) 



in terms of the extended Hamiltonian in (TT§j) . Evidently by f)20p and (1211) 
the above conditions i) and ii) are satisfied. Hence matrix G in ( 1201) is an 
example of a liner factorization of —K in ffl9|) . Other such factorizations of 
K include the non-symmetric matrices, 



G 



E-V bp 
ap 1 



for a ■ b = I /2m. However, it appears that K can be factored usually into 
Hermitian matrices. Thus for simplicity, though not necessary, it is assumed 
henceforth that any matrix factorization G of K is Hermitian as well as linear 
in the momenta. 

The fact that the determinant \G\ of G in (120]) is zero implies that G is 



a singular Hermitian matrix. But it is also well known in general that any 
matrix G is singular if and only if there exists some nonzero vector 6 such 
that 

G6 = 0. (22) 



Thus, in particular, the relation K 
relation, 

GO = 



E-V p/V2m~ 




'0i ' 


pj \Jlrn 1 







or that the two simultaneous equations, 

(E -V)9 1 + (p/V2~^j 6 2 
[p/V2^j 0i + e 2 



0, in ffl9|) holds if and only if the 

(23) 



= 0, 
0. 

0, (24) 



hold for some vector 9. 

Since K = 0, the most general solution of the pair of equations in ( 1241) is 
given by 



91 = A, 

9 2 = -(p/V2m^\, 



(25) 
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where A is an arbitrary function of x, t,p, E. Thus all vectors, 



9 



A 

Up/ \/2m) A 



(26) 



such that A is an arbitrary function, satisfy the null solution in (j22]) . 

It is demonstrated in the next section that the Fourier transform of 9 
in (1231) with respect to p is, in fact, closely related to the ^-function of the 
standard Schrodinger equation of a particle in a potential well. Thus the 
fact, that the factorization process of K always yields a singular matrix G 
and a nonzero solution 9 of G9 = 0, helps to supply an underlying basis for 
the Schrodinger -^-function of quantum mechanics. 

2 Fourier-Momentum Transform of GO = 
Using The Planck-de Broglie Laws 

In this section the Planck and de Broglie laws, 



imply that the relation, G9 = 0, in ( l22l) belongs to the domain of frequencies, 
uj and k. This suggests that the energy and momenta in ( 122]) be Fourier trans- 
formed to the spatial domain. Here in ( 1271) h = h/2ir, where h is Planck's 
constant, u is angular frequency in radians per second, p is the momentum 
vector of some particle and k is the wave-number vector in radians per unit 
distance associated with the particle. To perform such a Fourier transform it 
is convenient to assume that the generalized coordinates in ([3]) are Cartesian 
coordinates, e.g., q%, ?2, • ■ • ,q n are the components of the position vectors of 
the N particles, so that n = 3N. 

Now note that the set iV^ of all solutions 9_ of the relation, G9_ = 0, in 
(122]) is called the null space of the Hermitian matrix G. Also by (127]) any 
solutions 9_ in Nq can be expressed in the abbreviated form, 



E = huo and p = hk, 



(27) 



= 9(p,q) = 9(hk,q), 



(28) 



where 



9= (Qi,--- ,Qn,q n +i) 



(29) 
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and 



P = (Pi, ■ ■ • ,Pn,Pn+l), 

k = (h,--- ,k n ,k n+1 ) (30) 

are the generalized position, momentum and wave-number vectors with the 
(n + l)th coordinates given, respectively, by 

Qn+l = tj 

p n+ i = -E, and 

k n+ i = -u. (31) 

Since the matrix, 

G — G(p,q), (32) 

in (fT7|) is assumed to be the result of a Hermitian matrix factorization of K, 
all solutions 9_ of (|22|) automatically satisfy the constraint, K = 0, in (fTOj) . 
Hence vectors p and q in the function 9_ are independent vectors, e.g. see [1, 
Sees. 6.9 and 6.10]. This latter fact makes it possible to consider q to be a 
parameter vector of the system. 

Classically the equations of motion for vectors p and q are determined 
completely from the Hamiltonian by the Hamilton- Jacobi equations. How- 
ever, particles of a system, which satisfy the quantum rules in fT27|) . are 
wave-like in nature. Hence, the components of the generalized wave-number 
vector k in ( 1301) are the appropriate observables of the system. By ( 1281) any 
solution 9_ of (|22j) is a function of k and q as well as a function of p and 
q. Thus ^(p, g) can be considered to be some wave-number and frequency 
spectrum of the electro-mechanical system with q being a parameter vector. 

Explicitly this wave function is a (n + l)-dimensional Fourier transform 
of 9(p, q) of the form, 

±{q,l)= [ 9_{p,q)ey-UV p _, (33) 

where dV p is the volume element of the (n + l)-dimensional Cartesian space 
Rn+i- Here also q' is the position of the system, associated with this Fourier 
transform of function 9(p,q), and q is the position vector in 9(p,q) that is 
left unchanged by the mapping. 
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Now apply the same Fourier transform in (I33p to the relation in (j220 to 
obtain 

eksLGedVp = 0. (34) 

Rn + l 

Also remember that matrix G was chosen in the last section to be linear in 
the components of the momentum vector p. Hence to evaluate the above 
relation one needs integrals only of the form, 

/ Pk e^-6(p,q)dV p _ = ~-±{g,& (35) 



for k = 1, • • • , n + 1 where <j)(q, q') is defined in (|33l) . Finally a use of ([35 
and fl32l) transforms into 



G(l,q)£(q,l)=0, (36) 
where p' is the generalized momentum-operator vector, 

^ _h fd_ d \ 

- i \d<&' ' dq' n+ J ' 

To obtain a general Schrodinger-like wave equation which is a function of 
only a single position vector, one sets q' = q in ( 1361) . This yields 



G(£,q)0(q,£ 
or the vector wave equation, 



= 0, (37) 

<?'=<? 



G(p,q)±{q) = 0, (38) 
of the position vector g and the quantum-operator vector, 

t =U° » V (39) 



where 

±(q) = 0(q,q), (40) 

is a vector, Schrodinger amplitude function, of the n + 1 components of g. 
The general formula in ( 1381) for a wave equation is applied next to the matrix 
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factor G defined in (1201) of Example 1. 



Example 2 The use of (1201) in fl38l) yields the operator-matrix identity, 

= (41) 



E-V P/V2m 
p/ v2m 1 



^2 



or the pair of coupled operator equations, 



(p/V2m)ijj l + ^2 = 0, 



where 



P 



h d 



and Pi = —E 



(42) 
(43) 



ft 9 

i dx ' ' i dt 

Solving the second equation in (jl2"j) for tp 2 and substituting into the first 
equation yields immediately 



(p/2m + V)ifa. 
or by (jl3l) the wave equation 

9 2 



Eipi 



dx 2 



+ V )ipi = ih 



dt 



(44) 



(45) 



for the Schrodinger partial-differential equation of second order of a single 
one-dimensional particle in a force field in standard form. 

The above procedures, used in Examples 1 and 2, to factor the extended 
Hamiltonian and to find its associated Schrodinger equation, are applied in 
the next section to factor the extended Hamiltonian of an iV-particle system 
in an electromagnetic field. Formula (1380 is utilized then to find a general- 
ization of the Pauli equation from one particle to iV particles. 



3 Factorization of the Non-Relativistic Ex- 
tended Hamiltonian for N particles 

The relativistic Hamiltonian for N charged particles with magnetic potentials 
is known, e.g., see [6, Lec. 2], [7, Chapt. 12] and [2, Chapts. 7 and 8]. It is 
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given by 



N 



m k c 



k=l 



+ U 



E. 



(46) 



where U is the scalar potential energy of the N particles due to external and 
interacting electric fields as well as other forces acting on the system, and q k 
is the cahrge of the fc-th particle. Here as well A k is the vector potential of the 
magnetic field impinging on the fc-th particle, including externally imposed 
and mutually interacting magnetic fields of the N charged particled due to 
motion. 

The non-relativistic Hamiltonian of the system can be found directly by 
classical methods or by approximating the relativistic Hamiltonian as follows: 

N 




QkA k ) 2 + 



m k c 



+ U 



+ u. 



k=l 



Subtracting H from the total energy E yields the desired nonrelativistic 
extended Hamiltonian nullity relation, 



-K=E-U 



- 1 

V- 



k=l 



2mi 



Qk . 

p, - —A. 



0. 



(47) 



for the iV-particle system, where K is defined in ( fl~5l) . The higher order terms 
in the above expansion of f T46l) can be used to find relativistic corrections to 

dSD. 

It is demonstrated in the ensuing discussion that the determinant of the 
following 2 N (N + 1) x 2 N (N + 1) Hermitian symmetric matrix, 

D 1N A N ~\ 



G 



(E - U)I 



Pi- 



31*1 



Pi- 



/2mi 



'?JV A JV 



/2mi 



\/2m N 





\/2mjv 



(48) 
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is given by 



det(G) = K 2 



(49) 



where K is the extended Hamiltonian relation, and / and represent the 
unit and zero 2 N x 2 N matrices. In (148]) Pj and Aj for j = 1,2, ■ ■ ■ , N are 
defined respectively by the Hermitian matrices, 



p- = & r p and Aj = a r A p 



with 



(50) 



(51) 



[crij,CT2j,cr3j\ 

being the spin-vector operators for N independent particles. These operators 
are defined by the commutation relations 



[O'lj, &2k\ 

W2j,cr 3k ] 

[C3j, CTlk] 



idjkcrij, 
i5jk&2j- 



(52) 



with 5jk the Kronecker delta. They are traceless and idempotent: 



a 



a 



■2.1 



a 



3/ 



(53) 



For a single spin, these operators can be represented by the standard 2x2 
Pauli spin matrices: 



0"! 
0"3 



1 

1 

-i 

1 

1 
-1 



(54) 



By flSTil) . flST]) and (154"|) one would obtain more explicitly 



P 
A 



P3 Pi ~ W2 

Pi + ip2 ~P3 

A 3 A l - iA 2 
A 1 + iA 2 -A 3 



and 



(55) 
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for a one-particle system. For N particles, the spin operators become 2 N x 2 N 
matrices, which we can represent by 



cry 



hi- 1 



I 2 ® • • ■ 7 2 

N—j times 

® <7i ® (la)®"-'" 




0"! 



(56) 



and similarly for c 2 j and a^j, where (8> denotes the tensor (or Kronecker) 
product, J n is the nxn identity matrix, and ct 12i 3 are the 2x2 Pauli matrices 
from (1541 . By inserting these definitions into fl50|) . we get 2^ x 2^ matrices 
for Pj and A,-, which give us the 2 N x 2 N sub-block matrices in matrix G in 



<$E§ for j = 1, 2, ■ ■ ■ , N. 

To find the determinant |G| of G in (l4*9|) first simplify the notation of 
matrix G in (148]) by defining the following 2 W x 2^ sub-matrices of G, namely, 



al = (E-U)I and 
H,- 1 ^ * 



\j2rrij 



P., 



(57) 



for j = 1,2, 



iV, where a is the scalar function E — U. Then matrix G is 



re-expressible in the simpler sub-block form as follows: 



G 



al Hi • • ■ Hj\[ 

E x I ■■■ 

H N ••• / 

(1I2N B 

B In2 n 



(58) 



where I 2 n and I N2 n are 2 N x 2^ and (N2 N ) x (iV^) unit matrices. Also in 
(1581) 5 and are, respectively, the 2^ x (N2 N ) rectangular matrix, 



B = [H 1 ,H 2 



H 



(59) 



and its conjugate transpose. 

Note now that matrix G in (|58|) has the form of the general nxn matrix, 



M 



A B 
C D 



(60) 
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where A,B,C and D are, respectively, k x k, k x (n — k), (n — k) x k, and 
(n — k) x (n — k) matrices. It is shown in Appendix A that if the sub-block 
matrix D in fl60l) is nonsingular, then 

\M\ = \D\ ■ \A- BD^Cl (61) 



is the determinant of the n x n matrix M. Thus is fl60l) if one lets n = 
2 N (N + 1), k = 2 N and makes the correspondences, 

A «-> aI 2 N, B <-> B, 

C ^ B\ D^I N2 n, 

then by (161]) one obtains 

|G| = | a / 2 AT - BB^\. (62) 
But by (EHD and ([57]) one has 

AT 

BB^ = J2 H 1- ( 63 ) 

Also by (J57J the ifj in fl63|) are evaluated by the matrix squarings, 

H j = [hijVij + h 2j cr 2j + h 3j a 3j ] 2 

= hp 2 N, (64) 



where 

is the square of the magnitude of the vector — ^ — - — ^-^J > f° r J = 
1, 2, • • ■ , N. A substitution of ( 1641) into ()63l) yields finally the desired deter- 
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minant of G in (j55p . namely, 



\G\ 




N 

E-U-Y—fp.-^-A, 

j=i j 



K 2 = 0. 



(66) 



Evidently the result in (1661) for the Hermitian matrix G in (j4"5j) satisfies 
conditions i) and ii) in Sec. 1 for G to linearly factor the extended Hamil- 
tonian K in (14 7p . Thus by the arguments in Sec. 2 the quantum wave 
equation, associated with K in (j4"T|) and its matrix factor G in (JHJ), can be 
obtained by a use of the formula in (13T1) . To accomplish this note first by 
( 166]) that since the matrix G is singular, there exists a nonzero solution vector 

" ' 'Eat' Iv'" 'Sat' *) of the ec L uation 



(67) 



for the momentum vectors, p = (p l5 p 2 , • • • ,P N , —E) in (15U1) and the position 
vectors g = (g x , • • • , q N , t) in fljp. 
Next for the vector in (IBTl) let 



0-i 



1JV+1 



be the vector of the (iV + 1) 2 -component vectors, for k — 1, 2, • • • , iV + 1. 



The vectors ^ fc in (1681) are generalized spinor vectors or more simply, spinors. 
A substitution of the vectors in (1681) of spinors into (1671) . where G is given in 
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(I48p . yields the following system of linear equations, 



1 (fV-fVK + e, 



V2m/v V c 



iJV+l 



(69) 



in the iV+ 1 spinors 1 ,9 2 ,--- , #jv+i5 defined in (1681) . Thus in terms of these 
N + 1 spinor vectors the extended Hamiltonian relations, K — 0, in f|T6|) is 
valid if and only if the linear system in (169]) has non-zero solutions for spinors 
^15^2' ' ' ' i@-N+i f° r au momenta p , ■ • • ,p and total energy E. 

The null relation in flBTl) for matrix G in fH5|) transforms immediately by 
the procedures of the last section into the wave equation, 

G(p,q)±(q) = 0, (70) 

in fl38l). where 



E= 'Etv' - ^) ( 71 ) 

is the generalized momentum-operator vector with 

h ~ d 

p = -V fc and £ = %h— (72) 
- fc z dt ( J 

for = 1, 2, • • • , N being three-dimensional momentum operator vectors and 
the energy operator of quantum mechanics. Finally a use of (172|) in (150]) and 
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then a substitution into (1701) yields the "transformed" linear system, 



V2mi 

1 If, 11 



^ (A = 

1 (A,-^/w)i 1+ :0, 



-JV+l 



(73) 



of (TV + 1) coupled partial-differential equations for the (TV + 1) spinors, 
defined in ( |68|) . where 

-Pfc = g_ ■ p k = -a ■ V k and 

E = «| (74) 

for = 1, 2, • • • , TV, are momentum-operator matrices and the energy oper- 
ator of quantum mechanics. 

To convert the system in (1741) of linear, coupled, partial-differential 
equations for the N + 1 spinors ip l7 • • ■ ,if> N 1 into a single second-order, 
Schrodinger-like, equation in the single spinor if) first solve the fc-th equa- 
tion of the last N equations for ip to obtain 

i = -v^( A -7^)i (75) 

for k = 2, 3, ■ • • , iV + 1. Finally a substitution of these iV solutions into the 
first equation yields the second-order wave equation, 



(76) 

for the N non-relativistic particles of the spinor wave function 
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Next the energies are found for the interaction of the magnetic fields 
with the intrinsic magnetic spin moments of the N particles. The techniques 
used to accomplish this are adaptations of the Levy-Leblond approach first 
developed in [4] and in more detail in [10]. 

To find the interaction terms of the intrinsic magnetic moments of the N 
particles with the magnetic fields experienced by the particles one needs to 
extend the following well-known identity. For a single spin- 1/2 particle, 



AB = a - P + ia ■ (ax P) 

where a is the Pauli spin- vector operator in fl3T|) . a 
[Pi, fli-, Ps] are two 3- vectors of operators, and 



(77) 



[ail, «3] and (5 



A 



B 



a ■ a 



0-3 ax — ia 2 
«i + iat2 — CI3 



(78) 



The identity in (JTTJ) is easily seen to apply in the iV-particle well, by 

the definition (T5^i) of the kth spin operator a, where now we define A k = ak-a k 
and B k = a k ■ (3 , and the identity becomes A k B k = a k ■ (5 + ia k ■ (a k x /3 fc ). 
If we apply it to the kth squared matrix of the last iV terms of operators in 
the partial differential equation in fl76|) . we get 



c 

+ia k 



—At 



x [P 



Qk A 

v k - -Ak 

- q -A k 



(79) 



By (1721) p k is a differential operator which does not necessarily commute with 
a function. Thus the term in the brackets when applied to the wave function 
Tp , is simplified as follows: 



P k x lk 



f(^xp,)-f(P,x^) + |(^x^) 



x h)tx~ q i [(p h x A k ) - {A k x p k 

Qk 



--(lk X Ak)±y 
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Hence a substitution of this result into (1751) yields the identity, 



Qk 



Ai 



qk 



Qk 



.Qk 



—A, 



Hk . R 
^k ■ B k 

c 



(80) 



for k = 1, 2, • • • , iV since V& x A fc = B_ k is the magnetic field in the vicinity of 
the fc-th particle. This result is the same as that obtained in [10, Sec. 13.2] 
for N = 1. 

Finally a use of (1501 in (175]) and a replacement of ^ by ip yields the 
following general Schrodinger-like equation, 



N 1 

E-U-Y-^-ip-^A, 
^ 2m k \ Lk c - 1 

k=l K 



N 



k=l 



2mk,c 



'Ak-B k )U = 0, (81) 



or its equivalent in differential-operator form, 



ih 



0_ 

Of 



N 1 
u-V — 



k=l 



?v 4 - ^ 

Z C 



iV 



fc=i 



[ok-B k ) U = 0. 



(82) 

For N = 1 the wave equation in (152"j) agrees exactly with the Pauli equation 
with spin, e.g., see [10, Sees. 12.5 and 13.2], where 



qh 
2mc 



a-B = pi-B 



is the magnetic interaction energy of B_ with the intrinsic magnetic moment 
of the particle being 



qti 



-a 



qh 



S 



(83) 



— 2mc mc 

in terms of the usual spin operator, S_ = (|)<r. Thus for N particles the total 
spin magnetic-interaction energy operator is given by 



A' 



E ■ 

-'-'spin 



N 



k=l 



2m k c 



(84) 



k=l 



in terms of the intrinsic magnetic moments B_ k , defined in ( 183|) . Finally if B_ 
is a uniform magnetic field, eq. ( 1821) reduces to the Schrodinger equation in 
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[9, Sec. 113] for an atom of N electrons. (Note that here we are neglecting 
dipole interactions between the different spins.) 

In the next section the general Pauli-Schrodinger wave equation for N 
particles is applied to the simple Zeeman effect for the hydrogen atom in 
a weak magnetic field. The problem is treated as a two-body problem, an 
electron and a proton, both with spin |, rather than as a one-body problem 
of a single electron with a nucleus of a large mass, e.g., see [10, Ex. 12.3]. 



4 The Weak Magnetic-Field Zeeman Effect 
for Two Oppositely Charged Particles 



In this section the splitting of the spectral lines in a weak magnetic field of a 
hydrogen-like system is studied using the generalized Paul equation in flHTT) 
for N = 2 particles. This new treatment yields the simple Zeeman effect 
with a slightly modified Larmor frequency which depends on the difference 
of the masses of the two particles. The simple Zeeman effect given here is 
similar to the classical treatment, which uses the Pauli equation for N — 1, 
in that the relativistic spin-orbit interactions are neglected, e.g., see [10, Ex. 



In the Zeeman effect the magnetic field B_ is assumed to be homogeneous. 
Hence there is no loss in generality to let only the z component of B_ be 
nonzero, i.e., 



where B is a constant. Since the vector potential A is related to B_ by 
B = V x A, the quantity, 



is chosen to be the vector potential associated with B_ in (1851) . 

The generalized Pauli equation in flHTl) or fl82|) for N = 2 has the form, 



12.3]. 



B = (Q,0,B) 



(85) 




(86) 



1 




) 



2 



(87) 



Hip 



2m\ 
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where 



for {k = 1, 2) and 



A,= [~By k ,^Bx k ,0) . (SS) 



B X =B 2 = B_ (89) 

with B_ being the uniform magnetic field given in (1851) . Eq. (IBB]) is used in the 
ensuing discussion to study the simple Zeeman effect of a hydrogen-like atom 
of two particles, the negatively charged electron and the positively charged 
nucleus. Here in ( IBB1) one assumes that the charges and masses of these two 
particles are, respectively, 

qi = -e, <? 2 = +e (90) 

and mi, m 2 , where e = |e| is the magnitude of the charge of an electron and 
mi < m 2 . Finally the potential U in flHBl is the Coulomb potential given by 

U = - Z A (91) 
r 



where 



III -la I (92) 



is the magnitude of the relative position vector r, given by 

r = r j — r 2 (93) 

with r_ k for k = 1, 2 being the position vectors of the /c-th particle. 

The spin degrees of freedom for this system have a 4-dimensional Hilbert 
space. This can be expressed separately in terms of the spins of the two par- 
ticles; but one can switch to a different representation in which the magnetic 
moments are parallel or antiparallel. We will simplify the equations by as- 
suming that the magnetic moments of the two spins are oppositely directed, 
i.e., we restrict ourselves to single two-dimensional subspace. In this case, 
we can treat the system as if it had a single spin- 1/2 degree of freedom. We 
will denote the effective operator for this spin by a. (This is energetically 
favored by the spin-spin interaction, though we are otherwise neglecting that 
effect in this section.) 

Since the magnetic field B_ in fl85|) is assumed to be weak, terms of order 
iifc = Ak ' Ak i n dHH]) are neglected, and one obtains the approximate wave 
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equation, 

Hip ' ' 



— [v\ + —Ai ■p 1 )+—[p 2 + —A 2 ■ p 2 
m y- 1 c - l j fi2 \- 2 c - 2 , 

Ze 2 1 d 
+ U = ih—ip, (94) 



2M 

ehfl 1V A „ 
+ — fr-B 



r ) — dt— 



where 

// fe = m k /M (95) 
are the normalized particle masses for k = 1,2, and 

M = m 1 + m 2 (96) 

is the total mass of the two particles. 

The position vectors r k of the two particles are transformed next into the 
relative position vector r in (1931) and the center-of-mass position vector R by 
the following two relations: 

r = r 1 — r 2 and R = ii\r_ x + fi 2 r 2 , (97) 

where /ii and /i 2 are the normalized masses of the two particles, defined in 
fl95|) . This transformation of coordinates, using f!72|) and the rules of partial 
differentiation, yields readily 

P 1 =P_ + HiE. and p 2 = -ft + fi 2 P (98) 

as the transformation of the momentum operators p and p 2 onto the mo- 
mentum operators p and P. Here specifically 

p = -/iV r (99) 

is the relative-position momentum operator in terms of the classical reduced 
mass /i, where 

111 1/1 1\ , 
- = — + — = — - + - , 100 
H m 1 m 2 M ^/ii /i 2 / 

and V L is the gradient operator with respect to the components of r. Also 

P = -MV R (101) 
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is the center-of-mass momentum operator in terms of the total mass M in 
(pp|) and the gradient operator Vr with respect to the components of R. 
Also a substitution of the relatoins in fl88|) produces the vector potentials, 



(102) 



A = A n -[iiA r 



12 — £±R 

in terms of the new vector potentials, 

B 



A 

A 



2 
B 



R 



(A*2 -Hi)[-y,x,0] 
(ji2-{Jii)[-Y,X,0]. 



(103) 



Evidently A r and A R are functions of the relative and center-of-mass co- 
ordinates, respectively. Finally a substitution of the momentum operator 
relations in and the vector potentials in (11031) yields the Pauli-like wave 
equation 

— (?+- [(^ - ^)Ar + Ar] ■ P] + ( (I 2 + ~)Ar ■ P] 
/il/X 2 \ c J \ C J 



2M 



+ 



ehB 



1 1 

Hi 



a — 



± = fait* 



(104) 



for two particles in a weak magnetic field in the relative and center-of- 
mass coordinates. 

Now express equation (I104p in spinor form as follows: 



Htp 



' Hi 











>1 " 





H 2 




. ^2 _ 


= lh dl 


. ^2 _ 



(105) 



where 



Ho 



H' r + H' R + 



H' r + H' R - 



ehB 


" 1 


1 


2c 


mi 


m 2 


ehB 


' 1 


1 


2c 


mi 


m 2 



(106a) 



(106b) 
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Here also from (11041) one has 



H' = — 



V + -(A*2 - A*l)Ar 



2/i L- c 



e 
r 



2/i 



(107a) 



and 



i 



2M 
1 

2M L 



P 2 + ^(^.£ 



(107b) 



where 



A„ 



4 R 



1 



-.4 



R 



Ar + 

~ Hz — Hi 
+ A r = VRT{R,r) 



(108a) 
(108b) 



are the gauge transformations of the vector potentials, A r and A R = as- 
sociated with the apparent interactions energies of the linear and angular 
momenta of the system with 



a(r;R) 



(r ■ A 



fJL 2 - Hi 

R-A, 



(109a) 
(109b) 



being, using (I104p . the corresponding gauge functions. 

In order to separate the variables in the vector partial differential equa- 
tions in ffT05]) let 

i>k = fk{t)g' k (R;r)<P' k (r-,R) (110a) 
be the factorization of ipk, where 



9'k = 9k(R)e*P 



—it 
Tic 
—it 

he 



r(R-r) 
a(R;r) 



(110b) 
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for k = 1,2 with fk(t), gk(R), i^k{f) being, respectively, functions only of t, 
R, and r. Also the gauge functions r and a are defined in ffTTO and ffl~09bl) . 
A substitution of (IllOap into (I105[) yields 

9 ■ 

H k fkg' k (t>'k = in -^fk9' k <f)'k = Wk)g' k <f>'k 



so that one obtains the ratios, 



HkgWk 

g' k 4>' k 



ihfk 
fk 



En 



(111) 



where Ej- is constant for all t,R,r and k — 1,2. Here is identified with 
the total energy of the two-particle system in the magnetic field B_. From 
the second equality in fill II) one obtains 

hit) = e-*** 

for the time function fk{t). 

A use of (JTD7aj) and fll07bl) in the terms of equation fillip , when rear- 
ranged, yields the equations, 



Hi — Et ) g' x $ x 



H 2 - E T ) g' 2 cj)' 2 



m + H' R - 



ehB 
~2c 



' 1 


1 " 




rrii 


m 2 _ 


— E T j 


' 1 

mi 


1 " 
m 2 _ 


-Et) 



From this pair of equations one obtains the ratios, 



frl _|_ ehB_ ( j. i_ 



E> 



T 



H' R g 



jjl _ ehB_ I J. 1_ 

21 2c \ mi m.2 



E> 



T 



H' R g' 2 



92 



El, 



E, 



(112a) 
(112b) 



where El is again a constant, but in this case for all R, r and k = 1,2. El 
is associated with the translational energy of the system. 

The second equalities in f lll2al) and ( lll2bl) yield immediately by f ll07al) 
the result, 



2M 



c — 



g'k + E Lg' k = o, 



(113) 
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for the motion of the center-of-mass of the two-particle system. By (lllObj) 
and the gauge transformation in (1108bj) one finds the identity, 



9 k 



-V fl + -V r t 
i — c — 



exp 



-le 



he 



tg'k, 



-ie 



he 



which transforms (I113P into the Schrodinger-equation pair, 



hf_ 

2M 



V 2 R g k + E L g k = 0, (k = 1, 2), 



for the translational motion of the two-particle system. See [10, Sec. 9.1] 
for an excellent discussion of the gauge invariance of Schrodinger's equation, 
except for the usually unwanted phase factor. 

In a similar fashion the first equalities in (1112 all and flll2b[) produce the 
wave equation pair, 



eB 



1 



2c Vrrti 



1 

m 2 



0"3 



(114a) 



for k = 1,2, where 



1 

2^ 
1 

271 



P+ -(^2 - Vl)A r 
L- C 



e 
r 



P 2 + -(/X2 -Vl){A r -P) 
- C - - 



e 
r 



(114b) 



is the gauge-transformed wave equation in the relative coordinates r of the 
two-particle system. Here also E = E? — E^ is the total rotational energy of 
the system. 

Next observe by fl86|) that the term A r ■ p in (I114ap is given by 



A r ■ p = ih 



B 



d 8 
^ dx dy 



B 



lis) 



where L z is the z-component of the angular- momentum operator, e.g., see 
[8, Sec 4.8]. Thus a substitution of f 1 1 1 5 j) into flll4bj) and thence into flll4al) 
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Ho<P + ^-(~ —J (4 + ^ 3 )0 = #0, (116) 



yields the wave equation 

e 

— 2C \ffli 777.2/ 

where 

H = f/2/2 + e 2 /r (117) 

and e is the magnitude of the charge of an electron. 

Now define a two-particle Larmor frequency as follows: 

eB ( 1 1 \ eB , . 

where is the "increased" mass, given by 

1 1 1 

777. £ 777i 777.2 



(119) 



needed in the orbital magnetic- moment energy operator, eBL z /2crriL. In 
terms of the new Larmor frequency in (j!18p the wave equation in fl 1 1 6 [) is 
expressed finally by 

H o + uj l (L z + ha 3 )0 = Efa (120) 
or in spin components by the pair, 

Hot + WLiLg + hadfa = Efa, (121a) 

#00 + ^(4 + ^)02 = E(j) 2 . (121b) 

If the magnetic induction B is zero, i.e., B = 0, then 

01 = 02 = Ipnlm = RnlYl m {9, 0) (122) 

are the well-known standard eigensolutions of the pair of equations (1121 ah 
and (1121bj) . e.g., see [10, Sec. 9.2], where here by fllOOl) mass m is replaced 
by the reduced mass /1. Note also that the wave function ipnim in (11221) is an 
eigenfunction of L z , i.e., 

L z ^nlm = hmipnim, (123) 

e.g., see [10, pp. 67 and 245]. 
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Now let E n i m be the eigenvalues of the eigenfunction ip n i m of H in (11171) 
so that one has 

Holpnlm = E n i m 1p n i m . (124) 

Thus putting 



■nlm 







(125) 



into the left side of (T120D yields by (11231 and (11241) 







+ w L (L 2 + h) 



Ipnlm 





[£ n z m + u L h(m + 1)] 







Hence in (|T25l) is a spinor eigenfunction of the two-particle wave equa- 
tion in (11201) with eigenvalue, 



Similarly 



E 



(i) 

nlm 



E nlm + uj L h(m + V 



Enim + uj L h(m - 1) 



(126) 
(127) 



is the eigenvalue of the spinor eigenfunction, 



—nlm 







The figure in [10, pg. 248] illustrates this two-fold splitting of the energy 
levels of a hydrogen-like atom for n = 1, 2. 

The above application of the iV-particle Pauli equation in (1841) to the 
Zeeman effect for N = 2 evidently reduces tot he corresponding problem for 
N — 1 with a reduced mass /i, except for the modified Larmor frequency ujl, 
given in (11181) . The fact that the Larmor frequency u>l i n (H18p is expressed 
in terms of an increased mass in f 1 1 1 9 [) . instead of the usual electron mass, 
appears to have been noted previously for two special cases. 

In 1952 W. H. Lamb in [11], investigated by a use of the Dirac equation 
for a relativistic particle in a field the effect of nuclear motion ont he orbital 
energy of theZeeman effect of a hydrogen-like atom. He found in his eq. 
(153) that the effective g- value for orbital motion, taking the motion of the 
nucleus into account, is given by 



9l 



m 1 — m 2 
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Lamb's effective g-value is equivalent to the inverse in (11191) of m^, the "in- 
creased" mass. 

Finally it is known that in a weak magnetic field there is no Zeeman effect 
in positronium, e.g., see [12, pg. 214]. This agrees with the result obtained 
by letting mi = rri2 = m be the common mass of both the electron and 
positron, so that by (11001) and (11181) /i = m/2 and u L = 0. Thus, ignoring 
higher-order relativistic corrections, positronium has no weak-field Zeeman 
effect. 

The above observations agree with the new, more complete, Larmor fre- 
quency obtained in (11181) . This in turn suggests the validity of the generalized 
Pauli equation in (1ST]) and the power of the techniques developed in Sections 
1 and 2 to obtain this equation. 

5 Concluding Remarks 

The techniques developed herein for the general Pauli wave equation with 
spin generalize the results obtained by Levy-Le Blond in [4], but also the 
results in [8]. In [8], a similar, but special 4x4 matrix is singular if and 
only if a certain relativistic Hamiltonian relation is zero for a free relativistic 
particle. A procedure, similar to that developed in this paper in Section 

1, is used in [8] to demonstrate the existence of vectro solutions 0(p,q) of 
G9 = in the null space Nq of G, where p = [f- ,Pi,P2,P3] is the energy- 
momentum vector of special relativity. However, in [8] the properties of Nq 
are found in considerably more detail by an actual construction, using the 
special structure of the 4x4 matrix G. 

The quantization procedure, generalized in Section 2, is applied in [8] 
to the equation, G9 = 0, to obtain Dirac's equation for a free particle. 
The Fourier-momentum transforms of the functions in the null space Nq 
are shown in [8] to be the general solutions of Dirac's wave equation. It is 
believed that the more general methods developed here in Sections 1 and 

2, using the extended Hamiltonian, can be generalized to other problems of 
relativistic quantum physics. 

Related to this is the possibility that a higher-order Hermitian matrix 
G can be found which incorporates a finite number of relativistic correction 
terms of the series for H in (|46j) . With such a larger matrix, it might be possi- 
ble to use the procedures of Sections 1 and 2 to find relativistically corrected 
wave equations for iV particles with spin. However, the extra complexity, 
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needed to find these equations, requires that this relativistic extension to be 
accomplished elsewhere. 



Appendix A 



The formula in (1611) to evaluate the determinant of the nxn matrix M in ( T5TJD 
in block form results from the following factorization lemma for partitioned 
matrices. 

Lemma 1 Let A and D be, respectively, k x k and (n — k) x (n — k) 
square block matrices with D being non-singular. Then the identity, 



h BD- 1 

I n -k 



M 



A 



A B 
C D 



BD^C 
D 



h 
D C I n -k 



holds, where B and D are, respectively, kx(n — k) and (n-k)xk rectangular 
block matrices. 

Proof: This lemma is verified by a direct multiplication of the three 
matrix factors in the identity. The determinantal identity in (refeq:3.13) 
is obtained immediately by using the fact that the determinant of a prod- 
uct of matrices is the product of their determinants and Laplace's rules for 
evaluating determinants. 
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